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Abstract 

We consider inflaton decay to final state bosons with self-interactions of moderate strength. 
We find that such final state self-interactions qualitatively alter the reheat dynamics. In 
the case of narrow-band resonance decay, where a quantitative analysis is possible, we show 
that these final state interactions regulate the decay rate. The phenomenon of parametric 
amplification is then effectively suppressed, and does not drastically enhance the decay rate 
and reheat temperature. Detailed applications of our results to realistic classes of inflationary 
models will be considered elsewhere. 
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Inflationary cosmology provides solutions to the flatness, isotropy and stable relic 
problems of the standard hot big bang. In these models the universe experiences a period 
of superluminal expansion during which its energy density is dominated by the potential 
energy of a scalar field (inflaton). Inflation ends when the inflaton enters the oscillatory 
regime during which we have a matter-dominated FRW-universe, after which the inflaton 
decays to relativistic particles (reheating). Reheating represents the crucial transition from 
the epoch of scalar-field dominated dynamics to a hot FRW universe. After reheating, the 
universe becomes radiation-dominated and its evolution is just that of the standard hot big 
bang . 

In the standard picture of reheating |], @] the effective decay of the inflaton occurs when 

r ~ H, where T is the one particle decay rate, and the reheat temperature of the universe is 

i 

Tr ~ 0.1(r) 2 (from now on M p i = 1 and all dimensionful quantities are expressed in these 
units unless otherwise indicated). In this picture, which we will refer to as the linear regime, 
perturbation theory is taken to be valid and the occupation number of the final state bosons is 
assumed to be smaller than one (for fermions this is assured because of Pauli blocking). It has 
recently been recognized within different approaches M -|31| that this picture is incomplete, 



and nonlinear effects can change it essentially, leading by parametric amplification to an 
enhanced decay of the inflaton (for a recent review see, e.g. M). There are two different 
possible regimes of such a parametrically amplified decay. In the first the decay occurs 
over many oscillation times of the inflaton, but experiences parametric amplification with 
moderately large occupation numbers for the modes of the decay product field; this "narrow- 
band resonance" case is amenable to analytic calculation, and we will be able to analyze the 
modifications to parametric amplification due to final state-self interactions of the decay 
products quantitatively. In the second case of "broad-band resonance" the amplification of 
the decay is so strong that there is explosive decay of the inflaton field on a time-scale not 
hierarchically longer than the oscillation time, with large occupation numbers for the modes 
of the decay product field; this scenario is more difficult to analyze quantitatively, but the 
physical mechanisms which we discuss are of sufficient generality that we expect that they 
alter the decay dynamics in this case also. 
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The consequences of a parametrically amplified decay could be significant, and pose major 
difficulties in the construction of viable inflationary models. If the decay is very fast, the 
final state particle modes have, in general, very large occupation numbers and are far from 
thermal equilibrium; these fluctuations have certain effects similar to very high temperature 



thermal corrections [L2| , and after thermalization may themselves lead to high reheat 
temperatures, with thermal energy densities of order the inflationary energy density. While 
on the one hand GUT symmetry restoration due to these non-thermal configurations revives 



GUT scenarios for baryogenesis |23], [16], |17j , on the other hand it reintroduces the problem 
of heavy topological defects whose solution was one of the initial motivations for inflation. 

Furthermore, realistic models seeking to implement inflation need to be able to stabilize 
the required flat potential against radiative corrections. The only presently known method to 
achieve this in the presence of gauge, scalar potential, Yukawa, and gravitational interactions 
is by using (approximate) supersymmetry to enforce the appropriate non-renormalization. 
Parametrically amplified decay of the inflaton, and the resulting efficient reheat, poses a 
mortal danger to realistic, super symmetric, inflationary models. In supersymmetric theories 



the reheat temperature is constrained |]32j , [33], |3J] by the need to avoid thermal overproduction 
of gravitinos. This leads to an upper bound on the reheat temperature of order 10 8 — 10 9 GeV, 
which is orders of magnitude below what would be achieved by efficient reheat from amplified 
inflaton decay. One obvious way to avoid this disaster would be to insure that the inflaton 
is sufficiently weakly coupled to its bosonic decay products that the decay never experiences 
parametric amplification. This is typically the case for inflatons in a hidden sector which 
have only gravitational strength couplings to their (observable sector) decay products. On 
the other hand, in many models one introduces direct superpotential couplings of the inflaton 
to the chiral scalars into which it decays, and we wish to consider the nature of the reheat 
dynamics in such models, and whether they are ruled out by constraints on reheating. 

In particular, we will examine the effects that the final state self-interactions of the decay 
products have on the parametric amplification of the inflaton decay. Because of the large 
occupation numbers for the modes of the produced decay bosons, we expect that the presence 
of self-interactions of these bosons will result in large effective masses being induced for these 
modes. If the bosons are thermalized these may be interpreted as thermal plasma masses 
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from self-interaction; more generally they will occur as induced effective plasma mass terms 
in the mode equations for the decay field. As the mode occupation numbers increase, so do 
these induced masses, until they equal the mass of the inflaton, cutting off the decay. Decay 
resumes with the thermalization of the decay products, and their dilution and redshift by 
cosmic expansion, such that their induced masses dip below the inflaton mass; the system 
thus proceeds in a quasi-stationary process of decay and dilution such that the induced mass 
of the decay products is always of order the inflaton mass. This regulates the parametric 
amplification of the decay, preventing abrupt and efficient reheat. In our discussion we will 
analyze the effects of final state self-interactions for inflaton decays that would otherwise 
be in the regime of narrow-band parametric amplification. However, because the dominant 
effect is a kinematical cutoff of the decay, due to the self-induced plasma mass of the final 
state decay products arising from their strong self-interaction, we expect similar effects in 
the case of broad-band resonance. Indeed, recent studies of the broad-band decay regime 
|7|, [T(|, [L4[] indicate that the explosive decay from the sequence of higher resonance bands can 
only effectively produce particles whose mass does not exceed that of the inflaton by more 
than an order of magnitude. So in this case too there will be a kinematical cutoff due to 
the final-state self-interaction induced plasma mass of the decay products; although it will 
now be regulated to be not more than of order ten times the inflaton mass the qualitative 
effects should be otherwise similar. Applications of our present analysis to realistic classes 
of supersymmetric inflationary models will be considered elsewhere |3T| . 

As a basis for our subsequent arguments we will consider a chaotic inflation model with 
the following potential, whose features we take to resemble the generic features of scalar 
potentials which arise in supersymmetric theories: 

V-=imV + W>X 2 + /*W + <?Y (1) 

where for schematic simplicity the inflaton and the matter scalar x are taken to be real 
scalar fields and the self-coupling of the \ field is considered to be of moderate strength 
1CT 1 < g 2 < 1. In supersymmetric theories where the decay scalars are standard model 
chiral scalars which are gauge non-singlets the quartic potential terms in x arise as D-terms 
and the coupling is of gauge coupling strength g 2 . For inflaton-scalar couplings arising from 
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superpotentials in supersymmetric theories one has a = 2hm, and the cubic and quartic cou- 
plings of the inflaton to \ are related to each other. In general the superpotential couplings 
h may be, and in viable supersymmetric inflationary models are usually chosen to be, much 
smaller than gauge couplings h ^ g. The inflaton mass m must be bounded by m < 1CT 6 
in order to be consistent with COBE data on microwave background fluctuations []35 . 



First let us review the effects of parametric amplification on inflaton decay, ignoring self- 
interaction of the decay products. The nonlinear effects that lead to amplified decay act in 
two different regimes: 

[1] ^ 4> ^ ~ f° r the cubic coupling and ^ < 3 < for the quartic coupling. This 
is the narrow-band resonance case which can be analyzed perturbatively, and the dominant 
effect is the large occupation number for x's. This case has been considered in [0, |], pi], 
11, , where it is shown that parametric amplification occurs and there are narrow-band 
resonances for x production at k = — and k = m for the cubic and quartic couplings 
respectively. This is the case for which we will make quantitative estimates of the effect of 
inclusion of final state self-interaction of the decay products. 

[2] <p > — for the cubic coupling and </> > ? for the quartic coupling. Here perturbation 
theory is not valid and problem is highly nonlinear. In this case there is broad-band resonance 
for a large domain of momenta [[7], [TD| that leads to an explosive decay of the inflaton. As 
noted above, because of the essentially kinematic nature of the cutoff we expect the final- 
state self-interaction effects in this case to be qualitatively similar to those in the narrow 
band case. 

Note that if the cubic and quartic couplings of the inflaton to the decay scalar are of the 
form arising from a superpotential coupling [a ~ 2hm], then the conditions for parametric 
amplification are more general for the cubic coupling and it will dominate the decay. Another 
important observation is that although most inflatons may decay during a stage of amplified 
decay this does not lead to the decay of the entire energy density of the inflaton. In the 
case of narrow-band resonance the decay stops no later than the time when <p has the 
minimum value in the abovementioned range (^- for the cubic coupling and (7^) 3 for the 
quartic coupling) J7], ||, while in the case of broad-band resonance it stops at the time of 
the transition from broad-band to narrow-band resonance when the decay becomes out of 
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equilibrium [|7|, |H], After that, the remaining energy density of the inflaton is redshifted as 
R~ 3 where energy density of relativistic x' s is redshifted as R~ A . The decay of the inflaton 
will then be completed as in the usual picture, and if the energy density of the inflaton 
dominates at that time there will be significant dilution of relic densities from the first stage 



of reheat |p5j jig] . In supersymmetric theories this dilution is not, generally, sufficient by 
itself to solve the problem of overproduction of gravitinos. However, combined with the 
effect of final state self-interaction, which we consider below, it can successfully resolve the 
gravitino problem in many models. 

Now let us consider generally the changes to the parametric amplification reheat dynamics 
that arise from the self-interactions of the final state decay products. Since by assumption 
the inflaton decays to observable sector standard model (s)particles the final state bosons 
carry gauge quantum numbers. For these fields self-interactions with couplings as strong as 
the gauge coupling arise at tree-level from D-terms in supersymmetric models [], 

and at the one-loop level in the non-supersymmetric case. In addition, there are couplings 
of non-singlet scalars to the gauge fields as well as possible large superpotential couplings. 
The decay of the inflaton produces quantum fluctuations along the direction of final state 
particles in field space and drives them to large field values where the effect of self-coupling 
becomes important. These large field values induce a self-mass for the decay products that 
subsequently slows down the decay and tends to shut it off . 

During the oscillatory regime 4> = (po cos (mt) , and immediately after the end of inflation 
4>o ~ 10 _1 . Subsequently <po decreases with time because of decay and Hubble expansion. For 
4>o ^ x ^ ne i n fl a t° n decays predominantly via the cubic term. As was discussed in |7|, || |I0|| , 
in the range (^) 2 ^ 4>o ^ t 1 parametric amplification occurs and there is a narrow-band 
resonance for \ production at = -Sr. In this range for the inflaton field we can use the 
Feynman diagrams for one-particle decay, but the occupation number of final state particles 



1 There are, of course, directions in scalar field space which are both D-flat and F-flat in the supersymmetric 
standard model; in realistic no-scale supergravity models they will in general already have developed Planck 
scale vevs during inflation [Q. If the decay coupling of the inflaton happened to align along one of these 
directions then final state interactions of the type we consider would not affect this particular decay mode. 
More generally the decay products will themselves, in turn, have decay modes along these directions. We 
consider these issues elsewhere [Bf|. 
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is non-trivial and must be taken into account in the calcualtions. According to [J7L Bl 10 
parametric amplification will effectively convert most of the energy density of the inflaton into 
radiation once it is in the afore-mentioned range. Assuming rapid subsequent thermalization 
this leads to a reheat temperature Tr ~ mh~^ that is much higher than that of the usual 
picture ~ O.lmi/t for reasonable values of h. The quartic self-coupling of x induce 
a finite temperature correction to the mass-squared of Xi if they are thermalized, of order 
g 2 Tfi 2 ~ g 2 m 2 h~ l at this time which is much larger than m 2 (as we mentioned above the non- 
thermal corrections that exist before thermalization are even larger). However at an earlier 
time td when the thermal (or non-thermal f\ ) correction is of order ^ the one-particle decay 
becomes kinematically forbidden (note that the thermal correction to mass-squared of is 
of order h 2 Tji 2 which is normally smaller than m 2 for h a typical Yukawa type coupling). 

The Hubble expansion that subsequently occurs redshifts the correction to the mass- 
squared of x as R 2 an d this causes further decay. As long as td < H^ 1 these successive 
steps of expansion, decay and (perhaps) thermalization continue. Eventually the decay is 
not effective enough to compensate for expansion and there is a delay before the remaining 
energy density of the inflaton is converted into relativistic particles as in the usual picture, 
and the decay of the inflaton is completed. 

Now let us perform a detailed mode by mode analysis of the effect of quartic self-coupling 
of final state particles. Consider the potential 

V = ^m 2 4> 2 + 2hm<px 2 (2) 

where (f> and x are both real scalars. By mode expansion of x we derive the following equation 
for each mode 

R 

Xk + 3— Xk + (k 2 + 4/im0o cos mt)xk = (3) 

where a dot denotes differentiation with respect to time, and in this equation for the modes 
associated with comoving wavenumbers, we are using the physical wavenumber k, where 



2 Whether the thermal or non-thermal correction should be considered depends on how rapid the ther- 
malization rate is. Details will be discussed shortly. 
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Rk = kcomoving with R the scale factor . For the moment we will ignore the effect of Hubble 
expansion, since it generally occurs over a longer time scale than that of the effects we 
consider. We will consider issues of thermalization and Hubble expansion below [for other 
treatments of the effect of Hubble expansion on the parametric resonance see [0, H, |TT|1. 



By choosing z = yt, and in the absence of final state self-interactions we derive a Mathieu 
equation for the modes of the x field 

k 2 4hm(j)Q 

Xk +(j^ + -^rcos2z)x k = (4) 

where prime denotes differentiation with respect to z. In the case of narrow-band resonance 
(in which we perform our calculations) the Mathieu equation has resonance solutions in the 
first instability band (y) 2 — Ahmcpo < k 2 < (y) 2 + 4hm(pQ. Modes in this band grow expo- 
nentially in time, which one interprets as particle production. The (slow) Hubble expansion 
will eventually drive the modes out of the instability band, but they spend enough time there 
to reach a substantial occupation number |L8j . This time if, can be approximately calculated 



from 5k ~ kHtb where k = y and 8k = 8h<po is the width of the first instability band which 
gives t b ~ 

The quartic self-coupling g 2 x 4 will induce an effective mass-squared rn 2 e ff for the x field 
and the Mathieu equation with the addition of this self-coupling will be modified to become 

„ k 2 + m 2 eff Ahm(j)Q 

+( { rn )2 +-^COs2^) Xfc = (5) 

In a background of isotropically distributed x' s ove r a narrow band of momenta (which 
is the case in cj) — > XX decay before thermal distribution is achieved) m 2 e ff ~ g 2 -^ to the 
leading order, with n x the number density of x's and E x their energy [|7|, 0, 0, [U| [14]. 



After thermalization, in a thermal background of temperature T, we will have the standard 
result m 2 e ff ~ g 2 T 2 . 

At t — when <p starts oscillating m 2 e ff = and there is resonance in the band (y) 2 — 
AhrrKpQ < k 2 < (y) 2 + Ahm(j)Q. With particle production in this band rn 2 e ff increases 
and resonance in this band stops when m 2 e ff = 8hm(p after which the number density of 
particles produced in this band remains unchanged. There will be, however, resonance in 
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the band (tt) 2 — 12/im</>o < k 2 < (tt) 2 — 4/im</>o that also stops when m 2 e ff increases by a 
further amount of 8hm(pQ. This incremental change of rn 2 e ff and smooth transition from one 
resonance band to the next one continues until m 2 e // = (y) 2 ? after that there is no resonance 
solution for physical states, i.e. states with k 2 > 0. If we take E x ~ (y) for particles in all 
bands a change of 8hm(p Q in m 2 e ff correponds to an increase in number density of x' s by 
5n x ~ ^psi. For the first band (f f - Ahmcpo < k 2 < (f f + 4/im0 o and for hfa < m the 
occupation number /™ required to shift the resonance to the next band is calculated to be 

1 „ .m. 2 , , 4hd) m 2 „ An 2 , . 

on Y ~ of™ x 47t( — x ~ ^ — =>• fs. ~ — (6) 

X (2vr) 3 2 2 ; Y g 2 2 g 2 



It has been shown |L8] that in the small amplitude limit the analytic result of the decay 
rate in the n-th instability band can also be derived from a physical process, n particles with 
zero momentum that comprise the classical homogeneous inflaton field annihilating into two 
final state bosons. In particular, for the first instability band and for <C 1 we can use 
the one-particle decay rate from the standard perturbation theory Y = Narrow-band 
parametric amplification has been shown to be indeed an induced process |T^] which means 



f° r fh ^ 1 the enhancement of decay rate by the factor (l + /fc) must be taken into account in 
the particle calculation. The time ti which is needed to reach fm = 1 can be approximately 
calculated as f] 

r 9 2hm 2 d> 

h m , n — 12 — 4 

hmfo 2 ~ —^2- =>h~ yr-r (7) 



8tt ^ 7r/i0o 



and considering the enhancement factor the time <5£ which is needed to reach fm. = is 



4vr 2 

irh(j)o'" g 1 



St ~ 3^ In — (8) 



This is the time for the resonance in the band (tt) 2 — 4hm(po < k 2 < (y) 2 + Ahrrufio 
to stop. For the next bands k 2 is smaller and this means that less phase space volume 
is available for decay products or, equivalently, the occupation number for those bands is 



3 Note that e _rtl ~ 1 - Tti for Tti < 1, which is valid here. 
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larger. Therefore, more time will be needed for production of Sn x bands with 



<r 



smaller k 2 , but not greatly so, as the larger occupation numbers are obtained rapidly due 
the the large coherent final state enhancement. A reasonable lower estimate for the decay 
time td to effectively achieve m 2 e // ~ ^ is 



^~^r^~^iiT^ ln — ( 9 ) 



m , 4-7T 2 
In — 

3/im0 o " 8nh 2 (f) 2 g 2 

In order to have physically realistic estimates we take g 2 = 10 _1 in our calculations. This 
leads to 

8 m .„ . 

St ~TT ' 715X2 10 

So resonance at each band of width 8hrri(f) ends in a time ~ ^and in an approximate time 
of 4fe ^ o s decay effectively stops. We notice that t d < H^ 1 for <p Q > (^) 2 . 

In this analysis the effect of self-coupling of x on the solutions of the Mathieu equation was 
considered only to the first order, which is reasonable for the case of narrow-band resonance 
h(f)Q < m. The key assumption in our treatment of the Mathieu equation in the presence of 
the nonlinear term was adiabaticity, i.e. that we can use the instantaneous value of m 2 e ff 
which is also legitimate since it changes over a time 5t ~ ^ which is greater than m -1 , 
the period of oscillations of 0, again because we are in the narrow-band resonance regime 
h<p < m. 

Our results show major differences from the simple parametric amplification case. The 
effect of large occupation number for final state particles is not that dramatic here because 
there is a whole range of resonance bands instead of a single one and the effect of the self- 
interaction of the produced particles drives modes out of the resonance bands much faster 
than the Hubble expansion. Consequently the leading effect that influences the decay is the 
self-interaction of the decay products, which stops it very early. 

So far we have not considered thermalization of decay products and the Hubble expansion. 

The temperature of the thermal bath after thermalization of x's is calculated from Q 

4 In general the number density of particles is not the same before and after thermalization but the energy 
density is conserved. 
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— (Qb H — Qf)T 4 = p v = n \E Y ~ (11) 

where , 9f are the number of bosonic and fermionic degrees of freedom, respectively. We 
take the number of degrees of freedom to be the one in the minimal supersymmetric standatd 
model (gs = 9f = 74) which leads to a temperature T ~ 9^^f- F° r this temperature, 
however, the correction to mass-squared of x is 

2 

m 2 e// ~ g 2 T 2 ~ ^ (12) 

which is much less than (y) 2 for g 2 = 10 _1 . This is just what we expected because rn 2 e ff ~ 
g 2 ^- will be smaller after thermalization when the number density decreases and the mean 
energy of particles increases. Therefore if thermalization is effective, it will lower m 2 e ff 
significantly which leads to further decay and thermalization. This sequence of decay and 
thermalization stops when g 2 T 2 ~ ^ that is at a temperature T ~ m for g 2 = 10 _1 if the 
sequence is completed within a Hubble time. 

Considering thermalization of decay products and the Hubble expansion, there are dif- 
ferent possibilities depending on the relation among different time scales involved in the 
problem, t osc ~ mr\ St ~ ^, t d ~ t th ~ Q and £ H ~ V^bfc" The most 

important thing is that t osc is considerably smaller than all other time scales as long as we 
are in the narrow band regime 4>o < Therefore changes in m 2 e ff caused by decay or 
thermalization (that can be considered as changes in the background) are adiabatic and our 
analysis is in principle valid, irrespective of the relation among td, t t h, tn- Regarding these 
time scales there are different cases: 

[1]- 5t < t d < t t h tH- This occurs when the inequalities (j) > 5 x 10~ 3 ^ (to have 

U & t^), 0o ^ 3 x 10~ 5 (to have t t h & t H ) and f < (fio < f (to have t d < t H ) are satisfied 

which necessarily means ? < 10~ 2 . In this case the following sequence of events happens: 

decay of the inflaton to x's, end of the decay, and thermalization of decay products to a 

temperature T ~ m, all in a time scale shorter than the Hubble time. Hubble expansion then 

5 t t i l = Tth^ 1 with T t h ~ ct 2 ^ is for an out of equilibrium distribution of x's. It is easily seen that for a 
thermal bath with T ~ ^ (which is the highest temperature that can be achieved) m _1 < t t h also. 
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redshifts T but for O ^ (^) 2 (assuming td < t t h ^ tn all through the way) this sequence 
keeps T ~ m. For </> < (^) 2 decay is no longer effective to compensate for expansion and 
we must wait until a later time when O ^ h 2 and decay is completed as in the standard 
picture. This second stage of decay dilutes the gravitinos that are produced earlier during 
the first stage. 

[2]- St < tth <U< t H . This occurs for 0o < mm[5xl(T 3 f , 3xl(T 5 ] and (f f < O < f 
which is possible only for ? < 5 x 1CT 3 . This case is similar to case [1], only because of the 
faster thermalization the temperature is higher and closer to the maximum possible ~ m. 

[3\-St <t d <t H < t th . This occurs when 0o > Max[5xl0' 3 f , 3xl(T 5 ] and (f f < O < f 
that is consistent only for ^ > 3 x 1CT 5 . In this case decay stops without effective thermal- 
ization in a Hubble time, so the distribution of x's is out of equilibrium. Hubble expansion 
redshifts m 2 e ff as R~ 2 and further decay compensates for this change. A thermal distribu- 
tion of particles is not achieved, however, unless tth ^ tH Q ■ F° r 0o ^ (x) 2 ^ ne situation is 
as in case [1]. 

[4]- t th <St<t d < t H . This occurs when O < 10~ 3 f and (f f < O < f which is pos- 
sible only for ^ > 10 -3 . In this case the decay products thermalize almost instantaneously 
and there is thermal equilibrium from the very beginning of the decay. Actually both 5t and 
td are much greater than ~ and ~ 2h Z^ Q s in this case because rapid thermalization of 
decay products in the thermal bath keeps the occupation number at each resonance band 
below one. This means that temperature can be much lower than its maximum ~ m and 
the situation is very similar to the one in the standard picture of reheating. 

Cases [1] and [2] are the most dangerous regarding the problem of gravitino production. 

In these cases thermal equilibrium can in principle be achieved from the very beginning and 

can last until the time H^ 1 = ^J^m~ 3 h 2 . Even in these cases the gravitino overproduction 

is not that serious since at most we have a thermal bath with temperature T ~ m for a 

time t ~ y^m~ 3 /i 2 , much shorter than t ~ m~ 2 (assuming m <C (^) 2 ) which is the case 

for a radiation-dominated universe with temperature T ~ m. Case [3], on the other hand, 

6 In this case Hubble expansion moves particles from one band to another one. This slightly lowers the 
time St spent in those bands because now there is an initial number of particles in each band. This, however, 
is negligible since Sn x is also redshifted as R~ 3 . 
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is the most secure since the distribution of x's is out of equilibrium, a distribution with less 
mean energy per particle and higher number density compared with a thermal distribution. 
Depending on the parameters of the model m, h and strength of the self coupling g 2 one or 
all of these cases can happen for (^) 2 < 4> < — but the temperature is at most of order m 
as long as the inflaton is in this range, and is redshifted after that. Also the thermal and 

2 

non-thermal corrections to the mass-squared of x are always < 

To compare our results with that of parametric amplification without the effect of final 
state self-interaction let us consider the case h = lCT 6 ,m = 1CT 7 , where the standard picture 
predicts Tr ~ 0.1m5/i = 10 8 GeV. In the simple parametric amplification case almost all 
the energy density of the inflaton is converted into radiation once </>o ^ x = anc ^ l eac ^ s 
to a very high reheat temperature Tr ~ 10 15 — 10 16 GeV which, from the point of view of 
gravitino overproduction, is a disaster. According to our present analysis we are in the case 
[3] in the above all the way from 0o = t — 10 _1 to (f>Q = (^) 2 = 1CT 2 . This means that 
thermalization is not effective for 1CT 2 < 0o < 1CT 1 and occurs much later when decay is 
no longer effective, so the temperature during the first stage of decay is actually lower than 
~ 10 12 GeV. This results in a gravitino number density after the first epoch of decay which 
is too large by a factor of at most 10 12 . The inflaton decay is then completed via a second 
stage when 0o ^ ^ — 10~ 14 . By this time the temperature of the thermal bath and the 
momentum of the relativistic particles that might have been produced during the first stage 
are redshifted by a factor ( ^pr ) 1 = 10~ 8 . The second stage will determine the effective 
reheat temperature to be ~ 10 s GeV, and releases a large amount of entropy that dilutes 

12 3 

the gravitinos produced during the first stage of decay by a factor of (^s-) = 10 12 which is 
now sufficient dilution factor. 

As we noted previously, since the effect which we consider is a kinematical cutoff due to 
the large self-induced plasma mass of the strongly self-interacting final state decay products, 
we expect to also have similar effects in the broad-band case. Studies of the broad band 
resonance case 0, [10], [HJ indicate efficient production of the decay products only for masses 
up to about an order of magnitude larger than the inflaton mass, so that when the broad-band 
resonance has built up a sufficient density of the decay products ( typically in a non-thermal 
"preheat" distribution) that their self-induced plasma mass exceeds this range, the decay will 
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be suppressed. As above the decay will subsequently proceed as thermalization and Hubble 
expansion reduce the number densities and the induced self-mass for the decay products dips 
into the accessible range, resulting in a regular, quasi steady-state transfer of energy into 
the decay products. 

Finally, we contrast the effects considered herein with those considered by Khlebnikov 



and Tkachev fl3fl , who studied the semi-classical non-linear effects of the inflaton decay 
coupling in a massless X(p 4 model. Note that the effects they study have a different and 
independent origin from those considered here. The effects which we study are specifically 
due to final state self-interactions of the decay products which are different from, and larger 
(gauge strength) than, the inflaton decay coupling. If these final-state self-couplings are 
present, then we saw above that their effects acted to regulate the parametric amplification 
of inflaton decay. 

Returning to the narrow-band CclSC, clS treated above, it is useful to ask if there are vi- 
able models where we are in the narrow-band regime from the beginning of oscillations, and 
for which the analysis presented here provides quantitative, as well as qualitative guidance. 
We note that in a simple chaotic inflation model with potential V(<p) = \m 2 (f) 2 inflation 
ends when V" (0) ~ 247rV(0) which happens for <fi < 10 _1 . In the case of primordial super- 
symmetric inflation or new inflationary models the amplitude of post-inflation oscillations 
around the global minimium is generically substantially smaller than the Planck scale. De- 
pending on the parameters then, some of these models may satisfy the inequalities necessary 
to be in the narrow-band resonance regime, or even for the nonlinear effects to be negligible. 
Interestingly enough, it seems that in some viable supersymmetric models this is indeed the 
case, and the analysis we have here undertaken is quantitatively valid. We will return to 
these issues elsewhere |H[ . 

In conclusion, we have seen that the self-interaction of final state bosons of moderate 
strength, that arises very naturally in supersymmetric models, has an important impact 
on the decay of observable sector inflatons, besides producing a rapid thermalization rate. 
As a step towards improved understanding of the reheating process we have considered 
a simple schematic model representing generic features of supersymmetric theories, with 
such a final state self-coupling, and have shown that in the case of narrow-band resonance 
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the outcome is qualitatively different from that of simple parametric amplification. Here 
inflaton decay occurs during two stages: one stage that consists of successive steps of partial 
decay, thermalization and expansion at early times which ends relatively soon, and a second 
stage as in the standard picture that completes the decay. In the first stage because of 
the quasi-adiabaticity we were able to show that the temperature and the amplitude of 
quantum fluctuations of final state particles are at most of the order of the mass of the 
inflaton (approximately), and the temperature is several orders of magnitude below the 
naive predictions of parametric amplification. The second stage of decay then determines 
the final reheat temperature and releases a substantial amount of entropy; is of particular 
importance in order to dilute the previously produced gravitinos in realistic supersymmetric 
models. 
Note Added 

After this paper was submitted a preprint appeared from Prokopec and Roos ||hep-ph /9610400 



who do lattice simulations of inflaton decay. In the case where the decay product field has 
final state self-interactions of moderate strength their simulations exhibit regulation of the 
broad-band resonance decay, in agreement with the general analytic arguments we have 
presented above for the broad-band case. 
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